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A System for Measuring 

the Thermal Resistance 

of Semiconductor Diodes 

A fast, automatic system for accurately measuring 
junction-to-case thermal resistance of semiconductor diodes. 

By Norman R. Galassi and Bernard S. Siegal 




Fig. 1. This complete thermal resistance measurement 
system speeds evaluation of diodes. In production, it al- 
lows fast checking for packaging faults. The diode under 
test is placed in the heat sink located in the center of the 
rack second from the bottom. 



HERMAL RESISTANCE OF SEMICONDUCTOR DIODES is 

1 one of the more important parameters specified on 
a diode data sheet. Knowing the guaranteed thermal 
resistance of a particular type of diode, a designer can 
determine if this diode will survive its own internally 
generated heat when used in the equipment he designs. 

The problem in making a determination, however, is 
that many published specifications are not guaranteed 
and often inaccurate, even if unintentionally. Inaccuracy 
arises because of a lack of standardization of measure- 
ment of thermal resistance, and because inaccuracies 
creep in due to rudimentary methods of measurement. 

Attempts have been made to standardize the measure- 
ment of thermal resistance. The Department of Defense 
and the Joint Electron Device Engineering Council 
(JEDEC) of the Electronic Industries Association (EIA) 
have published material toward this end. One such tech- 
nique, called 'Method 4081 of MIL-STD-750A; gives 
only an outline of a thermal resistance test. Another 
publication, the 'Technical Report on the Measurement 
of Thermal Resistance' of JEDEC describes several 
measurement methods and explains the sources of error. 

However, details to implement these methods are not 
fully described. Different manufacturers are using differ- 
ent techniques which lead to differences in their pub- 
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lished specifications. Therefore, as a safety factor, de- 
signers often reduce the rating from published specifi- 
cations by 50 to 90%. Diodes with small thermal time 
constants (e.g., some step-recovery diodes) are the big- 
gest problem since they heat and cool so rapidly that 
accurate measurements are difficult. To resolve these 
problems, Hewlett-Packard Associates has developed a 
system shown in Fig. 1, whose accuracy can be demon- 
strated. In addition, the system speeds measurement time 
from the 3 to 10 minutes now required to about 10 to 
20 seconds. 

For most circuit applications, the total thermal resist- 
ance that a device presents to its surroundings consists 
of several separate thermal resistances as expressed in 
Equation l. 1 



he + @cs + 



(1) 



where 

0j a is the total thermal resistance of the heat path from 

the diode junction to its surroundings, 
9j C is the thermal resistance of the heat path from the 

diode junction to the diode case, 
8 CS is the thermal resistance of the heat path from the 



diode case to the heat sink (mica washer, etc.), and 
6 sa is the thermal resistance of the heat path from the 

heat sink to the surroundings (ambient). 
The total thermal resistance (0 ja ) is the quantity that 
the circuit designer needs to know. Generally the semi- 
conductor manufacturer has direct control only over the 
junction-to-case thermal resistance (#,,), and this is the 
parameter usually specified on the data sheet. 

The thermal resistance of any diode or heat sink is 
defined as the ratio of the steady state temperature differ- 
ence between the particular points of interest on the 
device to the power causing the temperature rise. 2 This 
is expressed by: 

T T 

La 

where 

6j C is the thermal resistance of the particular heat path 

under consideration (junction to case), 
T c is the temperature of the cooler end of the path 

(case), 
Tj is the temperature of the hotter end of the path 

(junction), 
P a is the power applied to cause temperature rise. 




Fig. 2. The large thermal mass of the diode holder keeps the 
case temperature substantially constant during the measurement. 
The diode, right, is inserted into the small threaded hole. 




Fig. 3. A typical plot of junction temperature as a func- 
tion of time. Heating power is applied at T» and removed 
at T„. The slope of several thermal time constants, t u t«, 
t 3 , are shown. These may vary from about 100 /is to 10 to 
20 seconds. 




Fig. 4. Typical plots of temperature coefficient K (forward 
voltage at 1 mA) for two types of step recovery diodes. 
Curves 1 and 2 are for diodes of the same type from dif- 
ferent lots and show that K is essentially the same, but 
forward voltage can differ. Curve 3 is for a different type 
of diode, and both K and forward voltage differ from 
Curves 1 and 2. 



Measuring Thermal Resistance 

In order to determine 6 ic for diodes, it is clear from 
Equation 2 that three quantities must be measured: 

T c , the case temperature, 

Tj, the junction temperature, 

P a , the power applied to the diode under test in order 
to raise the junction temperature. 

To measure the junction temperature, it is practical 
to monitor a temperature-sensitive parameter. In semi- 
conductor diodes, the forward voltage drop at a small 
fixed forward current (0.01 to 5.0 mA) is an excellent 
parameter to monitor since it is fairly linear with tem- 
perature and is repeatable. 

If the diode under test is mounted in a fairly good 
'infinite' heat sink (9 ct -f- sa small), Fig. 2, the junction 
temperature before the heating power is applied will be 
the same as the case temperature. The case temperature 
changes little during junction heating. 

It is extremely important to remember that the steady 
state temperature difference is used in calculating thermal 
resistance. Most diodes exhibit several separately identi- 
fiable thermal time constants, Fig. 3. 3 This requires that 
the power be applied for a time much longer than the 
largest time constant of the diode so that the steady state 
condition is reached before the junction temperature is 
measured. In addition, it is necessary to interrupt the 
heating power to measure forward voltage at the forward 
current used to determine junction temperature. This 
must be done quickly after interrupting the heating 
power; the time interval must be small with respect to 
the shortest thermal time constant of the diode under test. 

To measure junction temperature by using the forward 
voltage of the diode under test we must first measure the 
temperature coefficient of forward voltage of the diode 
under test. The temperature coefficient is: 



K = 



Ti 



V fm2 - V 



/ml 



I fm = constant 



(3) 



where K is the temperature coefficient (°C/mV) 
7\ is the lower test temperature (usually 25 °C), 
T 2 is the higher test temperature (greater than the 
highest that could be reached during the measure- 
ment of 9 jc ), 
V fml is the forward voltage at I fm , T u 
V fm2 is the forward voltage at I fm , T 2 , and 
Ifm is the fixed forward measurement current. 
While there is a fairly large (±10%) variation in the 
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Fig. 5. Temperature coefficient measurement system (a) is used for characterizing new 

diode types or spot checking during production if desired. The system (b) (pictured in 

Fig. 1) determines thermal resistance in about 10 to 20 seconds per diode, which enables 

fast measurement of large quantities of diodes. 



absolute value of the forward voltage at a fixed tempera- 
ture and forward current for diodes from different lots 
(batches), the temperature coefficient is relatively con- 
stant (±3%) from lot to lot for diodes of the same type 
(number), Fig. 4. This allows us to measure the tempera- 
ture coefficient of a sample of diodes of a given type and 
use the average value of the temperature coefficient in 
the thermal resistance calculations for that type diode. 
To heat a semiconductor diode, a known fixed forward 
current is applied to the diode and the resultant forward 
voltage measured. The dc power thus applied is: 

P a = I Pl X Vfl, (4) 

where 

P a is the heating power applied to the diode, 

I fh is the heating current applied to the diode, 

Vf h is the forward voltage of the diode at I fh . 

The measurement of the thermal resistance of diodes 

then requires three steps. 

1. Determine the temperature coefficient of forward 
voltage (K) at the measurement current (I fm ) over 
the specified temperature range (7\ and T 2 speci- 
fied). 

2. Measure the required parameters — V fmj , V fmc . 
and Vf h at a predetermined I fm and //». 



where 

Vf m j represents the forward voltage corresponding to 

the temperature of the junction heated and 
Vfmc represents the forward voltage corresponding to 

the temperature of the case. 
3. Calculate the thermal resistance by: 

/i _ I y fmj V fmc)f*- 



v fh X / 



fh 



(5) 



Thermal Resistance Measurement System 

As shown, the measurement of diode thermal resist- 
ance consists of measuring several different quantities, 
then calculating the thermal resistance of the diode under 
test. This can be a very time consuming task, taking as 
long as three to ten minutes per diode. If an adequate 
characterization of the temperature coefficient is first 
made, Fig. 5(a), the system described in this article can 
calculate the thermal resistance of an average diode in 
between 10 and 20 seconds. This makes the thermal 
resistance calculation system quite useful for evaluating 
packaging techniques, package uniformity, and product 
quality in R & D, production, and Quality Assurance. 
This system has the following pertinent specifications: 
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Fig. 6. The thermal resistance tester portion of the measurement system contains the diode 

under test. Programming and timing of the power supply, left, is controlled by the thermal 

resistance calculator and sequencer section. 
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Heating current (/ /;i ) capability: 0.05 to 5.0 amps 
(Power capability to 40 watts — dependent on the 
diode 

Measurement current (//,„): 1.0 mA 
Measurement duty cycle: 1.0% 
Accuracy: +8, —10% 
The HPA thermal resistance measurement system con- 
sists of three main subsystems: 

■ The Temperature Coefficient Measurement System 

■ The Thermal Resistance Tester (measurement 
circuits) 

■ The Thermal Resistance Calculator 

Temperature Coefficient Measurement System 

In practice samples of diodes can be measured on a 
periodic basis to insure that the proper value of K (as 
defined in Eq. 3) is used in measuring the thermal 
resistance. 

The temperature coefficient measurement system used 
for this measurement, Fig. 5(b), consists of an oven, a 
block of large thermal mass (on which the diodes are 
mounted), a constant current source for supplying I fm , 
a temperature measurement system, and a forward volt- 
age readout, such as an HP 3440/3443 DVM. The 
diodes to be measured are mounted in the block and 
placed into the oven. The oven temperature is adjusted 
to the specified T u and the forward voltage V fml is meas- 



ured on all of the diodes at I fm and T t . The oven tem- 
perature is then adjusted to the specified T 2 , and the 
forward voltage V fm2 is measured on all of the diodes at 
If m and T 2 . The temperature coefficients are then deter- 
mined using Eq. 3, and the average value is calculated 
for a number of sample diodes. 

Since we have found the temperature coefficient K 
to be fairly constant over a wide range of temperatures, 
it might be more practical to use constant temperature 
water baths at 0°C and 100°C to make temperature co- 
efficient measurements on large quantities of diodes. 

Thermal Resistance Tester 

The thermal resistance tester, Fig. 6, performs the 
measurement of V lmc , V fmj , and V fh . V fmc and V fmi are 
measured during a 10 p,s interval when the heating cur- 
ent is interrupted and only I fm flows through the diode 
under test. V !h is measured during the time the heating 
current is applied to the diode under test. These functions 
are performed by four main sections: 

1. Shunt switch. The shunt switch, Fig. 7(a), removes 
the heating current I fh from the diode under test 
and the blocking diodes isolate it so only the 1.0 
mA measurement current I fm flows. Thus, V Jmc and 
Vf m j may be measured during the 10 /xs interval. 
The shunt switch turns off 7 /7i for a 10 /xs interval 
at a 1.0 kHz rate (1.0% duty cycle). 



2. Sample gates. The sample gates, Fig. 7(b), are used 
to measure V fmc , V fmj , and V fh . V !mc and V lmj are 
measured by Sampling Gate No. 1, but at different 
points in time (determined by the Thermal Resist- 
ance Calculator). Sampling Gate No. 2 is used to 
measure Vp,i 

3. Timing circuits. The timing circuits synchronize the 
sampling gates with the shunt switch so they switch 
on and off at the proper times, Fig. 8. 

4. Heating current source. The heating current source 
supplies the heating current I {ll to the diode under 
test. The heating current source consists of an HP 
Model 6266 Power Supply programmed to supply 
a predetermined current at a time determined by 
the Thermal Resistance Calculator. 

Thermal Resistance Calculator 

The Thermal Resistance Calculator, Fig. 9, acts as the 
main control unit for the Thermal Resistance Tester, 
current source, and DVM in addition to performing the 
operations required to solve Eq. 5. It consists of seven 
subsections: 



1 . Two sample and hold circuits. One circuit samples 
Vfmc (before the current source is turned on) and 
stores the value in analog form. The second sample 
and hold circuit samples V Smj (after variable delay 
interval after the current source is turned on so 
steady state T , is reached) and stores the value in 
analog form. 

2. The subtraction circuit. Performs the operation 
Vfmj —V fmc on the outputs from the sample and 
hold circuits. 

3. The power calculator. Performs the operation V {h 
times the programmed I fm to obtain P a . The value 
of P a is converted into digital form by the HP 
Model 3440/3443. 

4. The division circuit. Performs the operation: 



fmj 



fmc 



K 



Pa 



This is accomplished using a multiplying digital- 
to-analog converter and a comparator to divide 
(Vfmj — Vfmc) in analog form by P a in digital form. 
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Fig. 7. The shunt switch (a) of the thermal resistance tester section removes the heating 
current In from the diode under test. The blocking diodes isolate the diode under test from 
the shunt switch so only the measurement current flows through it. Transistors Q, through 
Qi should be capable of switching 1-2 amperes each in less than 0.5 fis, and the diodes D, 
through D, should have a minority carrier lifetime of less than 0.2 fis at 5 amperes. 

During the measurement interval, the sampling capacitors in the sampling gates (b) are 
connected to the diode under test through the FET's, which have a low 'on' resistance. The 
capacitors retain the value of the voltage sampled since the FET's have a high 'off resist- 
ance at other than the measurement interval, and the high-impedance buffer amplifiers 
isolate the measurement circuitry from the capacitors. 
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Fig. 8. Timing diagram of the 
thermal resistance measurement. 
Sampling Gate #1 samples Vt™ 
and it can be delayed as long as 
6 us from the turn-on time of the 
shunt switch. This delay permits 
diodes with a long minority car- 
rier lifetime adequate time to 
reach \ !m from V fh . Sampling 
Gate #2 samples V fh . 
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5. The timing circuits. These circuits time all other 
circuits in the Thermal Resistance Calculator as 
well as the power supply and the DVM, Table I. 

6. The programming circuits. These circuits program 
the power calculator as well as the power supply 
and the DVM. 

7. The digital readout. This consists of an HP Model 
3440/3443 DVM. The DVM first reads the value 
of the power input to the diode under test, converts 
it into digital form, and returns this data to the 
Thermal Resistance Calculator. The DVM then 



TABLE I Time Sequence of Operation of 
Thermal Resistance Calculator 



Time 


Action 


To 


Initiate Test 


T, 


Sample and hold V fmc 


T 2 


Apply lf h to the diode under test 


Tto 


Time delay* 


Ta 


Sample and hold V fmi ; calculate 
V fmi — Vfmc; calculate power 


T< 


Read power on DVM 


T 5 


Calculate e » V *"' ~ V,ml! ; 
Turn off U. 


T. 


Read 9 on DVM 



* The time delay between T 2 and T 3 can be set between 1 and 
30 seconds to allow the diode under test to reach thermal equi- 
librium before the thermal resistance is calculated. 



reads the thermal resistance with a scale factor of 
1.0 mV= 1.0°C/watt. 

Errors 

The main sources of measurement error are: 

■ Errors due to variations in the temperature coeffi- 
cient K (estimated to be ± 3 % ) . 

■ Errors due to thermal time constants caused by: 

a. Errors in measuring Tj before steady state is 
reached. These can be reduced by adjusting the 
variable time delay TD. 

b. Errors due to junction cooling during the meas- 
urement interval. These could be as large as 
— 5% for a device with very short constants 
(100 /ts) and close to zero for a device with long 
time constants. 

■ Errors due to inaccuracies in the Thermal Resist- 
ance Calculator circuits which could be ±2% for 
temperature rises of at least 60°C at the junction 
during testing. 

■ Errors due to T c (case temperature) drift during the 
measurement of Tj which could range from 0% 
to +3% for a well designed heat sink. 



Adding these gives a probable overall accuracy of -f-8, 
-10%. 
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Fig. 9. Thermal resistance cal- 
culator and sequencer programs 
and times the rest of the system. 
The relay switches the DVM in- 
put between P, and 9 circuitry. 
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Fig. 1. Pushbutton bank on 
the HP Model 5 105 A 
Frequency Synthesizer permits 
selection of high-quality 
signals up to 500 MHz 
in steps as fine as 
0.1 Hz increments. 



Digital Frequency Synthesizer Covering 0.1 MH; 



By Alexander Tykulsky 



Frequency synthesizers have become highly versa- 
tile instruments since their introduction as signal 
sources for traditional HF communications. They are 
applied as essential parts of sophisticated systems in such 
areas as spectroscopy, advanced radar, space communi- 
cation and automated testing. 

For some time, there has been a need for high quality 
synthesized microwave signals capable of being switched 
rapidly. Suitable synthesizers, until recently have barely 
reached the VHF range. The use of auxiliary multipliers 
has been the usual expedient for bridging the gap, but 
this is costly and does not lead to optimum performance. 

A new frequency synthesizer, the HP Model 5 105 A, 
(Fig. 1) represents a step toward better frequency cover- 
age of the active spectrum. Frequency range is extended 
to 500 MHz, an increase of ten times the range of pre- 
vious HP synthesizers. Design techniques used produce 
these higher output frequencies with less noise than could 
be obtained by multiplying up from lower frequencies. 



The generation of signals in the microwave range is 
also simplified, since by starting at the higher output fre- 
quencies offered by the Model 5 105 A, less multiplication 
is required and phase noise is appreciably less. 

Synthesis System 

Operation of the synthesis system is accomplished in 
two stages. In the driver unit (HP Model 5110B), a 
family of fixed frequencies is derived from a frequency 
standard. In the variable frequency unit (HP Model 
5 105 A), the fixed inputs from the driver are selected as 
required, and further arithmetic operations are per- 
formed to produce the desired frequency. Fixed filters 
are used, and frequency changes are effected only by 
the operation of solid state switches. The method is 
'direct' allowing no possibility of an uncontrolled or in- 
correct output. A simplified diagram is shown in Fig. 2. 

The technique of generating a finely-resolved signal 
by means of an iterative sequence of mixer-divider units 
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Fig. 2. Output from the switched oscillator is processed, carried through two branches and difference 
mixed in the output mixer. Any frequency excursion occurring in the switched oscillator is cancelled out. 



>. to 500 MHz in 0.1 Hz Steps 



is borrowed from the HP Model 5100A. 1 The signal so 
generated, covering 30 to 3 1 MHz in decimal steps from 
0.01 Hz to 100 kHz, is added into the N (narrow band) 
branch of the synthesis system at Mixer 1 and is multi- 
plied by 10 in a doubler and quintupler. As a result the 
coverage contributed by the mixer-divider units as seen 
at the output is increased to 10 MHz and the steps range 
from 0. 1 Hz to 1 MHz. A search oscillator is included in 
the mixer-divider sequence to provide continuous tuning 
coverage or FM. 

A signal in the 30 to 39 MHz range, selectable in 
1 MHz steps, is added into the same branch at Mixer 2 
and contributes 10 MHz steps at the output. One strip- 
line oscillator of a set of five (covering 390 to 350 MHz 
in 10 MHz steps) is switched on as required. Simulta- 
neously the spectrum generator is switched to provide 
a complementary signal so that the difference (330 MHz) 
required to activate the N branch is always present. The 

1 Victor E. Van Duzer, 'A 0-50 Mc Frequency Synthesizer,' Hewlett-Packard Journal, 
Vol. 15, No. 9, May 1964. 



switched oscillator frequency is multiplied by 10 in the 
doubler and quintupler of the W (wide band) branch, 
resulting in the 100 MHz step of the synthesizer. 

Finally, a difference is taken in the Final Mixer and 
this difference signal is amplified and leveled. 

Drift Cancelling 

Non-synthesized signals (not quartz crystal derived) 
from the switched oscillators can be used because the 
synthesis system is in a drift cancelling arrangement. 
Cancellation of frequency errors is accomplished by 
channeling the frequencies through two paths as shown 
in Fig. 2. Since any UHF oscillator frequency excursion 
is carried through both paths, it cancels out in the output 
mixer. There is a difference in transmission time through 
the two paths, however, and the correlation is not perfect 
for rapid variations. The potential problem is avoided 
by using oscillators which have no appreciable high fre- 
quency modulation. 
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Fig. 3. Two of the switched oscillators are shown here 
(the two-layer board is disassembled). They were designed 
for low noise by placing the active element (a transistor) 
at a relatively low impedance point on the resonator and 
designing the resonators as folded transmission lines in 
a rigid structure to reduce microphonics. In addition, no 
high impedance point is brought out of the structure and 
oscillator switching is arranged so that resonator Q is not 
degraded. 




Fig. 4. Typical performance in switching from 399 to 
401 MHz. Lower trace is the switching signal. Sweep 
speed is 5 /j.s per division indicating switching speed is less 
than 20 jis. In photo, the signals are beat with 400 MHz. 



Switched Oscillators 

Exploitation of the drift cancelling scheme depends on 
the availability of easily switched, high quality oscillators 
(Fig. 3). 

Resonators are tapped transmission lines folded 
within the laminated boards. Tuning is accomplished by 
isolating portions of one ground plane. Output coupling 
results from the proximity of the resonators to the 50 CI 
output strip line. No high impedance point is brought 
outside of the ground planes. The resulting compact 
structure yields the necessary high spectral purity, low 
microphonic susceptibility, and good long-term stability. 

Output Circuit 

Output level is detected at an internal point directly 
on the output transmission line and compared to a refer- 
ence voltage. Control current is developed through a 
temperature compensated dc amplifier and applied to the 
PIN attenuator associated with the final mixer. 

This method of leveling, while permitting efficient use 
of the amplifier output capability, makes the output 
power dependent upon the VSWR introduced by the 
combined effect of the load, transmission line, and con- 
nectors. At the calibrated output level of dBm and with 
an ideal load, the flatness of output over the frequency 
range and for specified environmental conditions is ± 1 
dBm. Harmonic distortion is typically less than 4%. 

External manual control of the leveling reference pro- 
vides an output capability of approximately —6 to +6 
dBm. A metering circuit indicates normal operation of 
the leveling circuit. 

Spurious Outputs 

Non-harmonically related spurious signals are more 
than 70 dB down from the desired signal. Spurious sig- 
nals close in to the desired signal tend to occur in PM 
sideband pairs. Those more than a few MHz away tend 
to be single, additive signals. 

Switching Speed 

Programmed switching may be obtained by use of the 
three remote control connectors in the back panel of 
the variable frequency unit. Switching speed is nominally 
20 /xs, depending upon the frequencies involved. One 
example using the switched oscillators and spectrum 
generator is shown in Fig. 4. 

Phase Modulation 

Phase modulation is accomplished at Mixer 1 by con- 
trolling the tuning of the output circuit. The usable mod- 
ulation rate extends from to more than 1 MHz, with 
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the deviation typically down 1 dB at 1 MHz. Variation 
of deviation with output frequency is less than 2 dB. 
Harmonic distortion is a function of modulation depth 
— at 1 radian 5% is typical. Sensitivity is nominally 3 
radians per volt. 

Phase Noise 

Spectral plots of phase noise at 50 and 500 MHz are 
shown in Fig. 5. Also shown are the results to be ob- 
tained by multiplying the HP Model 5 100 A in noiseless 
multipliers (X 10). It is clear that while the HP Model 
5100A is to be preferred below 50 MHz, the HP Model 
5 105 A is preferable to multiplication (by 10) of the HP 
Model 5100A. The advantage comes from the smaller 
magnification of independent noise originating in the 
filters of the driver unit. This is made possible by the 
use of the drift canceling scheme. 

Solid State Modules 

The design is based upon a large variety of solid state 
circuits, ranging from diode switches to S-band quintu- 
ples. Modular construction is used throughout to achieve 
the required isolation between circuits, to facilitate exten- 
sive production testing, and to simplify servicing. 
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Fig. 5. Comparison of spectral distribution of phase noise 
sidebands of the 50 MHz Model 5 100 A and the 500 MHz 
Model 5 105 A. The lower frequency synthesizer provides 
improved noise performance at the lower frequencies, 
while at higher frequencies the Model 5105A shows sub- 
stantially better noise performance than can be obtained 
by multiplying the lower frequency by ten. 



SPECIFICATIONS 



HP Model 5105A Frequency Synthesizer* 

OUTPUT FREQUENCY: 0.1 lo 500 MHz. 

DIGITAL FREQUENCY SELECTION: 0.1 Hz through 100 MHz 
per step. Selection by front panel pushbutton or by remote 
switch closure. Any change in frequency may be accom- 
plished in 20 /is typically. 

OUTPUT VOLTAGE 

FIXED: dBm ±1 dBm into a 50-ohm resistive load. 
VARIABLE: -6 dBm to +6 dBm into a 50-ohm resistive 
load, typical at 20°C. 

SEARCH OSCILLATOR/FM: A search oscillator provides con- 
tinuous variable frequency selection with an incremental 
range of 1.0 Hz through 10 MHz. Manual, or external volt- 
age ( - 1 to - 1 1 volts) control with linearity of ±5%. 
Maximum sine wave rate (external) is 1000 Hz, 

PHASE MODULATION: (50L> rear panel input) 43 radians 
maximum deviations; dc — 1 MHz rate. Typical sensitivity: 
+ 3 radians per volt. Linearity depends upon deviation 
amplitude. 

SIGNAL-to-PHASE NOISE RATIO: Measured in a 30 kHz band 
centered on the signal (excluding a 1 Hz band centered on 
the signal) is greater than: 



Output Frequency- 


-MHz 


1 


50 


100 


500 


Ratio — dB 


48 


48 


48 


40 



SIGNAL-to-AM NOISE RATIO: (Above 100 kHz): Greater than 
74 dB in a 30 kHz band. 



RMS FRACTIONAL FREQUENCY DEVIATION (with a 30 kHz 
noise bandwidth):** 



Averaging Time 


Output Frequency 


1 MHz 


50 MHz 100 MHz 500 MHz 


10 miliiseconds 
1 second 


1x10' 
2x 10 ■* 


2 x 10 ' 1 x T0-* 6x 10- |D 
4x 10"" 2 x 10" 1 x 10 M 



SPURIOUS SIGNALS: Non-harmonicaily related signals are 
at least 70 dB below the selected frequency. 

HARMONIC SIGNALS: 25 dB below the selected frequency 
(applicable to fixed output when terminated in 50 ohms). 

WEIGHT: Net, 75 lbs. (34 kg). Shipping, 133 lbs. (60 kg). 

PRICE: $9,150.00 

* Requires 5110B Driver. 

** Note: When trie 5110B Driver utilizes an external frequency stand- 
ard, this will 3ffect the stability and spectral purity of the 
output. Performance data stated above are based on internal fre- 
quency standard or indicate synthesizer contribution to over-all 
performance with external standard. 



HP Model 5110B Synthesizer Driver 

OUTPUT FREQUENCIES: Provides 22 fixed frequencies for 
the 5105A Frequency Synthesizer: 3.0 through 3.9 MHz in 
0.1 MHz steps (50 mV + 1. -3 dB) 30 through 39 MHz in 
1MHz steps, 24 MHz, and 20 MHz (100 mV ±1.5 dB), 50- 
ohm system. 1 MHz buffered output (1 volt ±1.5 dB into a 
50-ohm resistive load) available at rear panel connector, 

INTERNAL FREQUENCY STANDARD 

TYPE: 1 MHz Quartz Oscillator. 

AGING RATE: Less than ±3 parts in 10* per 24 hours. 
STABILITY: As a function of ambient temperature: 
±2 x 10" 10 per °C from 0°C to + 55°C. 



As a function of line voltage: ±5 x 10" 11 for a ±10% 
change in line voltage (rated at 115 or 230 volts rms 
line voltage). 

RMS FRACTIONAL FREQUENCY DEVIATION: 

Average Time 1 MHz Output Frequency 

10 milliseconds 6 x t0 -w 

1 second 1 x 10" 11 

SIGNAL-to-PHASE NOISE RATIO:* >85 dB. 

SIGNAL-to-AM NOISE RATIO:* >80 dB. 

HARMONIC SIGNALS: >40 dB below the output (with 
proper termination). 

* In a 30 kHz band centered on the carrier, excluding a 1 Hz band 
centered on the carrier. 

PHASE-LOCKING CAPABILITY: A voltage control feature 
allows 5 parts in 10 8 frequency control for —5 to +5 
volts applied externally. 

EXTERNAL FREQUENCY STANDARD 

INPUT REQUIREMENTS: 1 MHz or 5 MHz. 0.2 V rms min- 
imum, 5 V maximum across 500 ohms. Stability and 
spectral purity of 5105A Frequency Synthesizer will be 
partially determined by the characteristics of the external 
standard if used. 

GENERAL (Applies to 5105A/5110B) 

OPERATING TEMP. RANGE: to + 55°C. 

WEIGHT: Net, 54 lbs. (24 kg). Shipping, 67 lbs. (30 kg). 

POWER: 115 or 230 volts ±10%, 50 to 400 cycles, 35 watts 
each unit (independent supplies). 

PRICE: $4,350.00 

MANUFACTURING DIVISION: HP FREQUENCY & TIME 
1501 Page Mill Road 
Palo Alto, California 94304 
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Phase Noise in Frequency Synthesizers 



A FREQUENCY SYNTHESIZER OPERATES On the signal 
provided by a frequency standard to produce a de- 
sired frequency which reflects the accuracy and stability 
of the standard. In effect, the frequency of the standard 
is multiplied by a selectable synthesis factor — a rational 
number usually based on the decimal system. The ease 
and speed of selecting the synthesis factor, the resolu- 
tion and range of available output frequencies, and the 
spectral purity of the output signal are usually the critical 
characteristics of the instrument. 

Circuits performing the necessary arithmetic opera- 
tions introduce random amplitude and phase variations. 
Amplitude modulation tends to be reduced by subsequent 
non-linear, compressive circuits; phase modulation how- 



ever, remains, and is usually enhanced. Synthesizer noise, 
then, is mainly random phase modulation. 

The total phase noise of a synthesizer includes the 
noise originating in the frequency standard. A useful 
model, Fig. 1(a), shows a simple synthesizer as the com- 
bination of an ideal multiplier and a phase noise modu- 
lator acting on the output of a frequency standard. This 
results in a total modulation made up of two parts which 
can be treated separately: The frequency standard noise 
multiplied by the selected synthesis factor, and the sys- 
tem noise 4> s (t,F). The system noise results from the 
combined effect of the arithmetic circuits (mixers, mul- 
tipliers, etc.). The importance of a noise source depends 
on its magnitude and the multiplication which follows. 



(a) 



(b) 



(c) 



FREQUENCY 

STANDARD 

Frequency = w 

Phase Noise = o (t) 



1 MHz 

FREQUENCY 

STANDARD 

Frequency = Fw, 

Phase Noise = o (t) 



FREQUENCY 

STANDARD 

Frequency = «> 

Phase Noise = 0„(t) 



SYNTHESIS 
FACTOR 



SYNTHESIS 
FACTOR 

Fa 
390-400 



SYNTHESIS 
FACTOR 

Fb 
390-350 



SYNTHESIS 
FACTOR 

F 
(Fa-Fb) 
0-50 



SYSTEM 

PHASE NOISE 

0s«) 



SYSTEM 

PHASE NOISE 

«A(t,F) 



SYSTEM 

PHASE NOISE 

■%(tF) 



SYSTEM 
PHASE NOISE 



-► Frequency = Fw 

Phase Noise = F0 o (t) +0 s (t,F) 



390-400 MHz 




Frequency = (F A — F B )w 

Phase Noise = (F A - F B )0 o (t) + « A (t,F) - B (t,F) 



390-350 MHz 



50 MHz 
^ Frequency = F<o 

Phase Noise = F* (t) + s (t) 

Note: s (t)>0 s (t.F) in Fig. 1 (a) 
for the same value of F. 



Fig. 1. Synthesizer noise models. Narrow band system (a) is a basic form of a typical 

synthesis system. The wideband system (b) is a combination of two narrow-band systems. 

The wideband system simplified (c) in which s (t) = 0A(t,F) — B (t,F) and is virtually 

independent of the synthesis factor F. The numbers given are examples only. 
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Fig. 2. Distribution of phase noise components. These 
curves also illustrate the sideband distribution if the ordi- 
nate is considered to be two times the single sideband to 
carrier ratio, and the abcissa represents offset from carrier. 



Wideband output ranges are most readily obtained 
by down-converting the outputs of two higher frequency 
sub-synthesizers* each of which operates over a prac- 
tical, small percentage bandwidth. The larger synthesis 
factors of each branch lead to greater system noise than 
would be expected from a non-heterodyne technique. 

An example of a wideband synthesizer is shown in 
Fig. 1(b). The total phase noise modulation, <j> T (t), to be 
expected is 

(F A - F B )<j> (t) + F A <j> A (t) - F n <t> B {t) 

where the first term is directly related to the frequency 
standard, and the last two terms, constitute systems noise. 
F A and F B are the selected sub-synthesizer factors, and 
F A — F n is the overall synthesis factor. For simplicity, 
the model of Fig. 1(b) can be represented as shown in 
Fig. 1(c) and the system phase noise can be written 

£ r (0 = F+o(?) + fa® 

Typical spectral distributions for the model of Fig. 
1(c) are shown in Fig. 2, where <£> (/), $ lS (/) and <b T (J) 
are the Fourier transforms of the corresponding time 
functions. 

Noise derived from the frequency standard becomes 
significant as F is increased, and then predominates at 
very low modulation rates. The best frequency standard 

* Sub-synthesizers in this context are the synthesis systems, separate but related, 
contained within a frequency synthesizer. 



available, consistent with the cost of the synthesizer, 
clearly should be used. Generally, a crystal filter is em- 
ployed as part of the standard to improve the spectral 
purity. 

System noise is virtually independent of F and is the 
major contributor at rapid modulation rates. Location 
of the corner, fc, (Fig. 2) is dependent upon the narrow- 
est filter bandwidth in the system. In systems where very 
fast switching is not required, the corner fc may be 
moved in by the use of a phase lock loop as a narrow 
band filter, but at the cost of some increase of close-in 
noise. 

The proper design of mixing circuits is a critical 
problem. Low signal levels lead to poor signal-to-noise 
ratios and phase noise. High signal levels may introduce 
spurious modulation at discrete frequencies — an equally 
undesirable result. Mixers applied at the end of a syn- 
thesis sequence introduce less noise since no multipli- 
cation follows. 

Drift cancelling offers another opportunity to reduce 
system noise. The high sub-synthesizer frequencies are 
obtained by adding coherent signals into each branch. 
The non-coherent noise is reduced (since its synthesis 
factor is smaller) while the coherent noise is canceled 
in the output. 

The sub-synthesizer factors for Fig. 3 are about one- 
sixth of those shown in Fig. 1. The system noise will 
also be about one-sixth — if the noise associated with 
the oscillator cancels. A differential group delay through 
the filters in paths A or B makes the cancellation im- 
perfect. The net effect of the oscillator phase noise can 
be seen at final mixer. Arriving through path A the phase 
is <f> D (t + t), where <f, D is oscillator noise, t is time and t 
is group delay due to the filters. Through path B the phase 
is <f> D (i). The mixer output shows the phase difference: 

<f, D (t + r) - 4> D (t) 

This equals r^ D (t) if t is sufficiently small. Writing $«(/) 
as the Fourier transform of <f> D (t), /2tt/$(/) is the spec- 
trum of <f> D (t), and j2irfT(fe D (f) is seen at the mixer output. 

It can be seen that if the synthesizer system noise is 
to be minimized $ (/) must be small, particularly for 
large values of /. 

This reasoning also explains the need for a sharp filter 
in the frequency standard output. In Fig. 1(b), assume 
a probable differential delay, t, between sub-synthesizers 
A and B. Then an additional noise modulation T F B $oO) 
appears in the output and this leads to the additional 
spectral component /tF b 27t/$ (/). For sufficiently large 
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FREQUENCY 
STANDARD 

Frequency = Fwq 
Phase Noise = o (t) 



SYNTHESIS 
FACTOR 

Fa 
60-70 



SYSTEM 

PHASE NOISE 

*A<t,F) 




A B 

(t = Group Delay) 



SYNTHESIS 
FACTOR 



SYSTEM 

PHASE NOISE 

*B(t,F) 




390-350 MHz 



SWITCHED 

OSCILLATOR 

Phase Noise = <P D 



390-350 MHz 



0-50 MHz 

Frequency = 

(F A - F B ).»o 

Phase Noise = 
(Fa - F B )0o (t) + 
*A(t,F) - 0B(t,F) 



Fig. 3. Noise model for a wideband synthesizer employing drift cancelling. System noise 

is less than in Fig. 1(b) since F A and F B are smaller leading to smaller values of A and 

<I>b. The numbers given are examples only. 



values of /, this component becomes overriding and must 
be reduced by filtering. 

Many different synthesis schemes can be devised to 
obtain a given set of numerical results, but the perform- 



ance of each is unique. Choice of a suitable design is a 
compromise involving consideration of discrete spurious 
outputs, switching speed, size, cost and other character- 
istics in addition to phase noise. ■ 
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Transform Methods for Linear Systems 



By Michael O'Flynn 

Professor of Electrical Engineering, San Jose State College 



This is a highly condensed collection of reference material 
on transform methods. Formulas for direct and inverse 
Fourier, Laplace, and Z-transforms are given in compact 
tabular form, with examples of how they are used to find 
system responses for several types of inputs. 

Transforms are powerful tools for solving linear time- 
invariant system problems, such as 



Given the i 



and given 

H(s), the system functio 
nputx(f)-*-B or 

/)((), the impulse respon 
of the system 



:tion 



find the output y(f) 



Depending upon the system and the excitation, the most 
convenient transforms to use might be Fourier transforms, 
one- or two-sided Laplace transforms, or one- or two-sided 
Z-transforms. These transforms are treated in many text- 
books (such as references 1 and 3 in Bibliography), and no 
attempt has been made to give an exhaustive treatment here. 
To the author's knowledge, however, there is nothing in the 
literature equivalent to Table I, which is a comprehensive 
collection of the defining formulas for all of these trans- 
forms and their inverses. The arrangement and notation of 
this table are designed to demonstrate the similarities and 
interrelationships between the different transforms. A fea- 
ture of the treatment of inverses is a logical progression 
from the fundamental inverse formulas to their equivalents 
in terms of complex variables (contour integrals, residues, 
Laurent series). 

Table II applies transform methods to the problem of 
finding the output of a linear system for five types of input 
signals. In each case the transform used is the most appro- 
priate one for the type of input. 

These tables can be used as a reference by anyone work- 
ing in the field of linear systems or — along with a textbook 
— by students of linear system theory. As an aid to the 
student, Table III lists the impedances of the basic network 
elements R, L, and C for the five classes of inputs. 

Classes of Inputs 

In the tables and discussion that follow, five classes of 
inputs are considered. They are: 

i) A damped sinusoidal input x(t) = Ae'"^ cos ( w „f + <j, ) 

oo 
ii) A periodic input x(t) = ]T A„ cos (n ao t + <f>„) 

n = Q 

iii) An aperiodic input which has a Fourier or a Laplace 
transform. 

iv) A sampled input x*(t) known only at specified points. 



v) A random input, where x(t) is a member of an ergodic 
random process for which the autocorrelation function 
R x (t) and its Fourier transform SJjw) are known. 

As an example of the use of the transforms, the aperiodic 
case (iii) will be discussed in detail. The sampled case (iv) 
is similar to the aperiodic case except that Z-transforms are 
used instead of Laplace transforms. The first two cases, i.e., 
periodic and damped sinusoidal inputs, are well known. 

For random inputs (case v) the problem can be restated 
as follows: Given the input autocorrelation function R.,(t) 
and the system function, find the output autocorrelation 
function ^(t). With this modification, the problem is similar 
to the aperiodic case. 

Inverses by Complex Variable Methods 

The problem of finding the output of a linear system 
whose input and system function are known may be re- 
garded as an exercise in finding inverse transforms. Inverses 
can be found either by looking them up in a table of trans- 
forms, or directly, by computing them using the inverse 
formulas given in Table I. The direct approach requires a 
knowledge of complex variables and contour integration, 
but it gives more insight into the physical meaning of the 
solutions.* For example, the inverse expressions contain 
full information about the natural frequencies of the system. 

Linear System Output for Aperiodic Input 

The output of a linear system for an aperiodic input 
x(t) is 



y(t) 






a+i oo 

H(s)X(s)e st ds. 



Here the Laplace transform has been used. H(s) is the 
system function and X(s) is the transform of the input x(t). 
If the one-sided Laplace transform is used, and if Jordan's 
Lemma is satisfied (usually this means simply that 
H(s)X(s)->0 as s^x) then 

y(t) = 0, t < 

y{i) = sum of the residues of the poles of H(s)X(s)e st , t > 0. 

The complete response y{t) consists of two parts. One is 
the forced response y f (t), in which only the input frequencies 
appear. [The input x(t) is often called the forcing function.] 
The other part of y(t) is the transient response y tr {t). 

* If a transform has an infinite number of poles (e.g., a discontinuous function, such 
as _T1_) it may be difficult or impossible to apply the complex-variable formulas. In 
such cases, time-domain analysis may be simpler than transform methods. 



Reprints of this article are available. Write 
Editor, Hewlett-Packard Journal, 1501 Page Mill 
Road, Palo Alto, California 94304. 
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Table I Definitions of Different 
Transforms with their Inverses 



Fourier Transform 



Laplace Transform (Two-sided) 



Laplace Transform (One-sided) 



Z-Transform (Two-sided) 



F(jo>) 



= 1 Kt)e~'">dt 

J -CD 



/ 



F(/'cd) is defined for m real and 
exists if 



r 

*s -or 



F(s) = I f(t)e- at dt 

S = a + fat 

F(s) exists if 



/ 



F(s) = / j{i)e~ st dt 

Jo 

S = a -\- jut 

F(s) exists if 



CO 

F(z) = £ «n)r" 



f(t) | d/ exists 



*-* -or 



/(f)e st I <A exists 



*' 



f(t)e st I eft exists 



or 

F(z) = F*(s) 



z = c 



This integral exists in some allow- 
able region of convergence 



This integral exists in some allow- 
able region of convergence a > <ri 



This exists in an allowable annulus 
of convergence ri < |z| < r 2 






which may be evaluated directly 
or from tables 



f(t) = ~ f FiMe**d(M 

1 f" 
= -^-r / F(s)e st ds 



For f < 

«0 = -^j(f)F(s) e °<ds 

" C, 

(Sum of residues of\ 
poles in right half J 
plane. / 

For t > 



KO = -±r(f)F(s)e°*ds 



Sum of residues of poles in 
left half plane. 









Cl < CT„ <(T- 



/W 



- 77 / J a.-iao 



F{s)e st ds 



CT > <Ti 



For / < 



M = -£j<f)ns)e ,, dt 



For r > 



/(0 



£0 



F(s)e l, 'ds 



C, 

(Sum of residues of\ 
poles with real parts ) 
to right of o- . / 



For t > 0, 

M = -^(f)F(s)e<d s 



= Sum of residues of poles 
with real parts to left of <r . 

Note: Residues are defined for 
counterclockwise contours. Con- 
tour d is clockwise. The integral 
around C, is, therefore, minus the 
sum of the residues of the enclosed 
poles. 



= Sum of residues of poles of 
F(s).[C encloses all poles of 
F(s).] 

Note: As a consequence of the 
definition of the one-sided trans- 
form an acceptable <j„ is automat- 
ically to the right of the poles of 
F(.s). 



Kn) = -^r(f)F(z)z"-'dz 



where C is a circle of radius r„ 
such that r, < r„ < r- 2 . 

For n S 

f(n) = sum of residues of F(z)z n ^ 
inside C. 

For n < 

f(n) = sum of residues of F(z)z"^ 
outside C. 

Note: the formula for f(n) is the 
formula for the coefficients of the 
Laurent series for F(z) (see Ap- 
pendix). 






Im(z) 



Region 
of 
Convergence 




-Re(z) 
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Z-Transform (One-sided) 



f*(t) = £ f(n)S(t - nT) 

71 = 

F(z) = £ Mz* 



F(z) = F*(.s)\ =eST 



This exists in an allowable region 
of convergence outside a circle 
centered at the origin. 



C 



where C encloses all the poles of 

f(n) = sum of the residues of 

F(z)z»- 1 

Note: the formula for /(h) is the 
formula for the coefficients of the 
Laurent series for F(z) (see Ap- 
pendix). 




Re(z) 



Table III 

Impedance Table 




Type of Input /(f) 


Transformed Input 


V(s) 
Impedance of Network Element: Z(s) = .. , 




R 


L 


c 


i) Damped Sinusoidal 
A e-".' cos (m t + </>„) 

Special Cases 

a) Sinusoidal, a„ = 

b) Exponential, 
oi = 0, </>„ = 

c) dc case, <t„ =: 0, 
„,„ = 0, f, = 


A e>-h 

A e»*« 
A 

A 

/ 


.phasor 
/notation 


R 
R 
R 
R 


Ls\ . . 
\S — — <r„ + lao 

Ls\ 
Ls\ 

\S — —o- 

.v = 

short 
circuit 


1 
Cs 

1 
Cs 

1 
Cs 

1 

C5 


s = —<To + M> 

■S = /wo 
,$ = (7 (} 

j = 

open circuit 


ii) Periodic 

a-- 
^A,, COS (na t -f tj>„) 


Same as case i) using superposition. 

Each harmonic has an impedance associated with it as in (i) (a). 


iii) Aperiodic /'(') 


Fourier Transform 
/(/») = ji\t) e-'-dt 

J -J3 

it integral exists 


R 


Ls\ 


1 

Cs 


s = /„, 


Laplace Transform 


R 


Ls 


1 
C.v 


iv) Sampled i*{ri) 


System functions for sampled data systems are usually computed 
first as ratios of Laplace transforms and then converted to 
Z-transforms. 


v) Random i(f) 

Autocorrelation 
Function = Ri(r) 


Power Spectral 
Density 

Si(j») = 

/£(T)r'"* 

•J LC 


\/R= 


Vir-vi . / i 

' s = /■» 
Note: S„(;») = | Z(/,„) |* S 4 (/<») 



Only the poles of X(s) contribute to the forced 
response. The poles of H(s) give the natural fre- 
quencies of the system and contribute to the 
transient response. If we denote the poles of X(s) 
as s t , s, 2 , . . . , s n and the poles of H(s) as S lt S 2 , 
. . . , S m , and assume zero initial conditions, the 
response y(t) is 

y(t) = y f (t) + y tr (t) 

n m 

i - l i-i 

where r x denotes the residue of pole s t in 
H(s)X(s)e !St and 

Ri denotes the residue of pole S t in 
H(s)X(s)e st . 

If the poles of H(s) and X(s) are all simple 
and distinct, and if the poles of H(s) do not 
cancel the zeros of X(s) or vice versa, then all 
residues will be of the form 

r, = r t «*i* and R t = RiA* 

See the Appendix or reference 2 for methods 
_ of computing residues. 



If the initial conditions are not zero the com- 
plete response can be obtained to within m arbi- 
trary constants. Again assuming simple, distinct 
poles and no pole-zero cancellation, y(t) will be 

n m 

*') = !>« +Z^ s ' f - 

The constants A f are determined from the initial 
conditions, which are usually given as values of 
y(t) and its derivatives at / = 0. 

For many forcing functions and for most 
systems, the forced response y f (t) is the steady- 
state response. Examples of aperiodic inputs for 
which this is true are step functions, ramp func- 
tions, and suddenly applied periodic functions. 

However, if the time constants of the natural 
frequencies — poles of H(s) — are greater than 
the time constants of the input frequencies — 
poles of X(s) — then the forced response will 
decay to zero faster than the transient response. 
In some systems, such as those made up of only 
L and C, and no R, the transient response will 
not die out at all. The transient response, there- 
fore, may actually be the total steady-state 
response. 
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Table II Input-Output Relationships for a 
Linear System. System Function = H(s) 



Type of Input 


Transform of Input 


Transform of Output 


Output 


Damped Sinusoidal 

x(t) = A e"»* cos (a„t + </>„) 

In phasor notation: 


X(s) = jnA [e'M(s — s a ) 
+ erf*, $(s — s *)] 

So = <T + /«0 


Y(s) = H(s)X(s) 


yft) = y [tf (,y )eM, e V + H(s„*)e-i^',"] 
= Re[H( Sl[ ) A eMeV] 


X T = A e**> 

It is understood that this means 

x(t) = Re[A e%e».<] 


Y T = H(s )X T 


y(/) = /?e[fl r ( J( ,)Z r e'.«] 


Periodic 

00 

x(t) = ^A n cos (n«i u ( -|- 0») 

71 = 

In phasor notation: 


Transform of the nth harmonic 
X„(ja) = Tr/4„[e J '*»8(m — nu ) 
+ e-'>,8( M + /!«,„)] 


By superposition 

Y(jm) = £ H(j a ) X„(j*) 

n = 


71=^0 


X. = A*». 


Y„ = HQnmo)X n 


XO = f;«e[ff(/«o J „)A'„e'»-„'] 

« -o 


Aperiodic 
x{t) 


X(j m ) — J x(t)e-'^dt 

•s -00 

if integral exists, 
or 

*W = / x(t)e~ st dt 


Y(j a ) = H(j*)X(j<*) 

or 

Y{s) = H{s)X{s) 


or 

i f*a+ioo 

y(t) = TT / [ff(fWj)]e"</j 


Sampled 

X*(t) = f) x(«)8 (< - nT) 

n=-00 


GO 

X(z) = X x(n)r* 

n--oo 


Y(z) = H(z)X(z) 
where 

00 

H(z) = £ h(n)z-" 

n--oo 

where 

h(n) = unit pulse response 


y(n) = ~-(pH(;)X(z)z''-<dz 
C 


Random 

R,(r) = Input Autocorrelation 
Function 


Input Power Spectral Density 
S/j*) = J R x (r)e-^d T 


S, ,(!„,) = | »(/„,) |« S^a,) 


Output Autocorrelation Function 

#»(r) = ~ f [\ H(j m ) \*SM»)] e)™do> 



Appendix. Complex- Variable Formulas 

1. (7) /(z)^z = if /(z) is analytic on and inside the 

„ closed contour C. 

2. (j) g(z)dz = 2-rt] ^2 [residues of poles of g(z) inside C]. 



3. If f(z) is analytic on C 1 and C, and in the region between 
them, then f(z) can be represented in that region by a 
Laurent Series. 



If g(z) has an nth-order pole at z then 

1 r d n 
residue of z 



M = Y l A n (z-z )» 

n--cc 

where A n = — - (7) 7 — r dz. 

2tt; jr (z-z„)" +1 




(n-1)! 



rfz" 



7 i (z - Zo)Mz) 



Z = Zn 



C is any curve between C, and C 2 and enclosing z . 
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